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ABSTRACT 


Vehicle suspension system plays a critical role in transferring static and dynamic loads produced by the vibration of vehicle body 
and wheels and absorbing shock vibration caused by road roughness. Low damping provides a better vehicle mass isolation and 
it give a ride comfort and hard damping provides vehicle stability with good road holding. The need to enhance conflicting odds 
between ride comfort and vehicle stability makes the design of the suspension a significant part particularly for off- road vehicles. 
Passive suspension can’t mitigate tradeoff between ride comfort and vehicle stability, because it presents a high frequency vibra- 
tion. It is in this line various types of suspension, like semi- active suspension, active suspension and intelligent suspension have 
been developed to reduce this compromise need. This paper aims to describe in details different types of vehicle suspensions, 
their characteristics, and their working principles mode. It illustrates in details the magneto rheological fluid (intelligent fluid) 
properties, compositions, mechanical model of Magneto-rheological damper like Bingham model and Bouc-Wen model. It also 
reviews semi active suspension control strategies based on Magneto rheological damper, like skyhook, ground hook, sliding 
mode, fuzzy logic and linear quadratic Gaussian. Simulations shows that a combination of more than two control strategies pro- 
vide a better vehicle comfort and vehicle stability at the same time. 
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1. Introduction 


The main key responsibilities suspension is to ensure ride comfort and 
road holding. It transfers different load and ground forces to the vehicle 
body. A tradition vehicle suspension constitutes by a spring coil spring 
or air spring and a damping element. The coefficients of damping ele- 
ment and spring are selected based on comfort, road holding, and han- 
dling characteristics[1]. Even if conventional suspensions may thrive a 
trade-off between ride comfort and road holding as though their spring 
and damping coefficients cannot be adjustably tuned depending on driv- 
ing efforts and road situation. It means that good ride comfort and road 
holding is gotten only under the designed circumstances. The studies and 
investigations on intelligent suspension with the development of micro- 
processor, sensors and actuators in 1980s have showed a great impact 
on avoiding the trade off by controlling the damping characteristics of 
suspensions. Based on the mode of generation of the desired control 
forces, intelligent suspension systems can be categorized into two main 
types: active and semi-active. The active suspension often requires sig- 
nificant amount of external energy to produce the desired control damp- 
ing forces. It is achieved by using pneumatic or hydraulic actuators . 
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1.1 Types of suspension 


In automotive dynamic era, suspension system is classified into 
primary and secondary, the primary suspension ensures the con- 
nection between vehicle body and wheels and secondary type rep- 
resent the seat. There are two basic types of elements in tradition 
suspension systems. These elements are springs and dampers. The 
spring in a vehicle’s suspension system ensure the mount of static 
weight of the vehicle and storage of energy applied to the vehicle 
body. The damper ensures the dissipation vibration energy and 
monitor different forces from the road transmitted to the vehicle. 
Basically the role and form of a suspension remain the same re- 
gardless of the category of suspension or automobile. Primary sus- 
pensions will be classified into passive, active and semi-active 
suspension. 


1.1.1 Passive suspension 


As it shown in the Figure! passive suspension is a type of suspen- 
sion which implies the control of inherent odds between mini- 
mum-frequency and maximum-frequency vibration mitigation. 
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This kind suspension system has of two important elements an en- 
ergy evacuating element (shock absorber), which is the damper, 
and an energy-storing element, which is the spring. As these two 
elements can’t give energy to the system, reason why this type of 
suspension systems is named passive. Passive suspension systems 
are exposed to different tradeoffs when they applied to a huge vi- 
bration frequencies[2]. Passive dampers figure out the conven- 
tional shock absorbers mostly applicable to many vehicles. They 
have the fixed damping characteristics, i.e. their damping charac- 
teristics cannot be adjusted while the car is being driven and thus, 
they are designed to achieve the most convenient trade-off be- 
tween ride, road holding, and stability[3]. 


Xs 


Sprung mass 


. 


Figure 1 Passive suspension 
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A high damping suspension will lead best vehicle road handling. 
However, it transmits much of the road forces to the vehicle body. 
Whenever the vehicle is riding at low speed on a harsh road or at 
high speed in a straight line, this will be perceived as a tough ride. 
The harsh ride is obnoxious for the drivers because it can damage 
the payloads or cause an accident. A soft damping suspension will 
lead a more comfortable drive, but may substantially mitigates the 
stability while the car turning, changing the lane, or exiting the 
ramp. Best design of a passive suspension can adjust some exten- 
sions on improving ride and road handling, but cannot remove this 
odd as indicated in Figure 2[28]. 
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Figure 2 Passive suspension trade off 


An ideal suspension system should provide good vibration isola- 
tion, i.e. minimum acceleration of the body mass, and ensure the 
maximal acceptable relative displacement between the vehicle 
body and different suspension components. 
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1.1.2 Active suspension 


The active suspension is a kind of suspension where the pas- 
sive shock absorber (damper) or both the passive damper and 
spring are replaced by an actuator force as it shown in the Figure 
3. Actuator force has the capacity to the spring and the damper to 
add and dissipate energy from the system, inversely to a passive 
damper, can only dissipates energy. In an active suspension, nec- 
essary displacement and velocity is applied independently by the 
actuator force. In active suspension the active force is provided by 
actuator, which is monitored and controlled by a control algorithm 
utilizing data from sensors mounted to the vehicle[4]. The active 
suspension has two important parties’ actuator and spring. The 
spring has a role supporting the static load of sprung mass and 
force actuator provides the necessary reactive force to mitigate or 
assimilate deflection caused by road abnormal. The actuator may 
be hydraulically, electromagnetically, pneumatically, or hybrid 
systems. These actuator systems are controlled and monitored by 
electric drives, the electric power of these drives is provided by 
the battery of vehicle or any other system in vehicle. Based to the 
high force ability, the easy design, advanced of technology, and 
commercial availability of the many parts, the hydraulic systems 
are commonly utilized in control of sprung systems of vehicle. 
The active suspension system equipped with electromagnetic ac- 
tuator is named electromagnetic active suspension system; it is 
constituted by spring and electromagnetic actuator fixed in paral- 
lel combination between vehicle body and unsprung mass. The 
electromagnetic actuator is controlled by embedded systems 
which supply the required electric power of the actuator. It gives 
active controlled force to mitigate road shocks in small time, elim- 
inates the roll and pitch vibration, and optimize road handling and 
ride comfort. When the actuator is operated and mounted mechan- 
ically it provides a high bandwidth active suspension and control 
body mass and unsprung mass. Low-bandwidth is provided to ac- 
tive suspension monitoring the vehicle body (sprung mass) when 
the active actuator operates and mounted mechanically in series 
with the spring and the damper. Generally, the range frequencies 
of the unsprung mass vary between 10-15 Hz, and the range of 
sprung mass frequency vary between 1-2 Hz. The desired fre- 
quency bandwidth determines the price of active suspension. 
When the desired bandwidth is limited the cost of control force 
actuator is minimize as cost of active suspension is minimized 
also[5]. With high damping force of active suspension, the best 
ride control is gotten. Active suspension operates by constantly 
sensing various vibrations in the road surface and provides those 
signals to the ECU, to activate actuator components. Actuators re- 
act upon the system to change its behavior, adjusting shock ab- 
sorber, spring rate and the like, to optimize ride comfort, drivabil- 
ity, and responsiveness. The ECU ensure the collections, analyzes 
and interpretations of data provided by sensors in little time like 
10 milliseconds. Nowadays a new kind of active suspension are 
being developed with linear electromagnetic motors monitored 
and controlled by microprocessors and complex mathematics 
equations[6]. 
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Figure 3 Active suspension 


1.1.3 Semi- active suspension 


Semi-active suspension systems have been investigated and stud- 
ied in the early 1970’s, the tradition spring element is remained in 
this kind of suspension, however the shock absorber damper is 
replaced with an adjustable and controllable damper as it is shown 
in Figure 4 below. 
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Figure 4. Semi-active suspension 
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Meanwhile an active suspension system requires an external en- 
ergy source to activate an actuator which is controls and monitor- 
ing the vehicle. The semi-active system utilizes just external 
power sensors and microcontroller to change the damping charac- 
teristics of shock absorber. The embedded microcontroller pro- 
vides the level of damping force based on a control strategy pro- 
gram to change the shock absorber characteristics to achieve that 
the necessary damping force. The semi-active suspension damp- 
ers adjust their damping force in real time according to a controller 
strategy which is often based on the system dynamics. By adjust- 
ing the damper’s fluid resistance movement, or damping coeffi- 
cient, changes the damper force. Based on its low cost and the 
capacity to change the damping coefficient independently of 
damper velocity, within limits, has increased a number of institu- 
tions companies and industries to investigate the ability of optimi- 
zation suspension performance by utilizing semi-active damper 
technology[7].The implementation of semi-active control law, is 
based on real time adjustment of damping force. Nowadays, many 
semi-active dampers are using hydraulic or electromagnetic sys- 
tems. The first type of semi-active utilizes mechanical valves 
which are activated by coil or stepper motor to control damper 
force in a hydraulic damper. The last kind of system uses control- 
lable fluids like electro-rheological or magneto-rheological to 
change their rheological effect in order to provide the required 
damping force. Even if the cost of semi-active suspensions re- 
mains high compare to the passive suspensions, they are very 
cheaper compare to active suspensions; as such, they are consid- 
ered as a more commercially viable option for many vehicles cat- 
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egories including commercial vehicle. Many researches and stud- 
ies have shown that semi-active presents the best performance of 
vehicle dynamic. 


2. Magneto rheological fluid 


Magneto-rheological (MR) materials are classified in smart mate- 
rials, they are presently being applied many applications. Smart 
materials are the materials which physical properties can be ad- 
justed or changes by external physical properties such as, temper- 
ature, pH, moisture, stress, electric or magnetic fields. Magneto- 
rheological material has various form of fluids; it can be a gel or 
even a solid material like elastomers. A magneto-rheological fluid 
mainly composed by micron-sized iron element that is immerged 
in carrier oil. Magneto rheological fluid presents the capacity of 
varying from free-moving liquid state into a semi-solid condition 
with limited fluid movement in fast response within some milli- 
seconds when applied to the external magnetic field. During the 
middle of the 20th century, researchers and scientists started to 
work with fluids able of experiment physical changes in a little 
time approximatively 10 milliseconds, for example, rheological 
fluids change viscosity or became semi-solids when are subjected 
to magnetic fields or electric currents. As it said above their vis- 
cosity varies in little time which is less than ten milliseconds. In 
four decades ago, magneto-rheological fluids were introduced in 
automotive applications. The main effort were emphasized on au- 
tonomous suspensions from the first intelligent approaches[8]. In 
1940s Jacob Rabinow [9] invented the fluid , which yield stress 
and viscosity changes when magnetic field is applied on it; since 
its invention different applications using MR fluid have been in- 
vestigated and designed such as dampers, clutches, engine mounts, 
haptic devices and building dampers, etc[10]. Physical character- 
istics changes of MR fluid are resulted from the chain-like struc- 
tures between paramagnetic MR particles in the low permeability 
solvent[11]. At the normal state, MR fluid presents the isotropic 
Newtonian behavior because the MR particles move freely as 
shown in Figure 5 (a). However, when the magnetic field is sub- 
jected to the MR fluid, MR fluid forms a chain structure as it is 
indicated in Figure 5 (b), MR fluids present the anisotropic Bing- 
ham behavior which result on resisting on flow and external shear 
stress. From this property, force or torque of application devices 
can be easily monitored by the density of magnetic field[12]. 
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Figure 5 Magneto-rheological fluid 


(a) Without Magnetic Field (b) With Magnetic Field 


2.1 Rheological properties of magneto-rheological fluid 


Rheology is the field which study and investigate the flow and 
changes of materials under exposed forces which is often meas- 
ured by a rheometer. The process of measurement and test of rhe- 
ological properties is applicable to all materials. Rheological 
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properties can be measured from bulk sample twist utilizing a me- 
chanical rheometer or on a micro-scale by using a micro capillary 
viscometer or an optical technique such as micro-rheology. Vari- 
ous mostly utilized materials and formulations show complex rhe- 
ological properties, whose viscosity and visco-elasticity can 
change according to the external conditions subjected, such as 
temperature, pH, moisture, stress, electric or magnetic fields. Rhe- 
ological properties is also measured in medical application where 
it given by internal changes of biopharmaceutical element such as 
protein concentration , stability or any kind interaction of between 
its properties[13]. Characteristics effects of magneto rheological 
fluid is based on increasing or decrease the viscosity of a magneto- 
rheological fluid and presenting yield stress based on a magnetic 
field applied on it as indicated by chain form of iron micro sized 
particles[14]. The MR effect is controlled and monitored by mag- 
netic field intensity and rheological characteristics of magneto- 
rheological fluid components. It means that rheological properties 
depend on particles polarity induced magnetic field applied on 
MR fluid. When there is a presence of a magnetic field, every mi- 
cro iron particle is changed into a dipole and makes a chain struc- 
ture with its consecutive particles which can resist failure for some 
ratio of heat and thus form a semi-solid structure. The interactions 
between the induced dipoles provoke the micro iron particles to 
be aligned along exposed field and make a form of column struc- 
ture. The chain structure will then prevent the fluid motion and 
consequently rise the viscosity of the suspension[15]. Among 
other factors, the rheology of MR suspensions depend on, particle 
shape distribution, quantity of particle, properties of the constitu- 
ent oil, additional additives, applied magnetic field, temperature 
[16]. In normal application, the off-state behavior of MR fluids 
relies on carrier fluid properties, additives, particle volume frac- 
tion, etc, as it shown in table 1, while the on-state performance 
depends on the solidity stage properties and the volume fraction 
of the solid stage. Rheology of MR fluids is explained in terms of 
pre-yield as well as post-yield conditions[17]. 


Ve = 0,T <T, 
T2To 


() 


Where, Ve indicates wrench and G the post-yield behavior of MR 
fluids provided by experimental test. 


Table 1. Table of rheological property of MR fluid 


Property Typical value 
Initial viscosity 0.2-0.3 [Pa.s] at 25°C 
Density 3-4[g/cm?] 
Magnetic field strength 150-250[kA/m] 
Yield point to 50-100[kPa] 
Reaction time 10 milliseconds 
Typical supply voltage 2-25V, 1-2A 
and current 


Work temperature -50 tol50 [°C] 


Characteristically, the diameter size of iron particles varies be- 
tween ranges of 3 to 5 microns meters. If the MR fluid is made 
with big size particles; it becomes difficult to stabilize the sus- 
pended particles. The cheaper iron particles widely available on 
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different market are carbonyl iron; they are usually inhibited to 
sizes bigger than | or 2 microns. Smaller size particles provide an 
easy suspension. However, the manufacture process of such par- 
ticles remains complex. Important smaller sized iron particles 
available on the market are usually oxides only, like pigments 
popular used in magnetic recording media[18]. In initial applica- 
tion where there is no presence of magnetic field MR fluids are 
estimated as Newtonian fluid[19]. Due to low voltage power sup- 
plies for MR fluids and relative temperature stability between — 
40°C and +150 °C, MR fluid present a great advantages more at- 
tractive smart materials than Electro Rheological fluids. MR flu- 
ids do not show yield stress, but show a variation of viscosity. 


2.2 Application of MR fluid 


Important use of MR fluids are in devices which require rapid con- 
tinuity and reversible variation in rheological characteristics. 
Magneto-rheological devices showed a significant advantage dur- 
ing the last decade, since MR fluid links mechanical parts indi- 
rectly with electronic embedded controller where the control strat- 
egies are incorporated on it, thus determine continuous mechani- 
cal characteristics of device. The use of magneto-rheological flu- 
ids is advanced of dampers, clutches, brake, control valves, vehi- 
cle power steering devices, artificial joints, alternators, sound 
propagation, chemicals sensing applications and other application. 
Especially shock absorber are the commonly device which utilizes 
MR fluid[20]. 

Valves: When the magneto-rheological fluid flows into the valve 
and magnetic field is subjected to the fluid provoke the increment 
fluid viscosity. This variation of viscosity generates resistance 
against flow of fluid in the valve. Thereby the pressure is raised 
and flow reduces the speed considerably or stops completely. 
Rosenfeld et al investigated structure and performance of mag- 
neto-rheological and electro-rheological valves and the result 
showed that performance of the MR fluid significantly depends 
on driving force and active volume of the fluid[22]. 

Polishing: Optical polishing, it is a good application of MR fluid. 
In working operation when a magnetic field is available, magnetic 
particles make chain-like structure and provoke abrasive particles 
to increase the polishing work. The ability of abrasive to keep 
power and stiffness of a MR fluid is greatly dependent on the size 
of induced particle. Type of abrasive and concentration are the 
other parameters which increase the quality of piece surface. 
Brakes: Magneto-rheological brake is an element which is mainly 
used in torque transmission. The torque of rotation MR fluids 
brake is changed quickly according to the magnetic field applied 
on it. Based on its simplicity design, ease control, the use of MR 
fluid remain a best choice for braking systems[23]. 


Dampers: Damper is an apparatus which is used to attenu- 
ate shock vibration, in addition, to mitigate vibration motion, and 
to minimize the kinetic force. A damper is called magneto-rheo- 
logical damper when it contains magneto-rheological fluids and 
monitored by with variation of magnetic field. They have been 
designed and developed to provide damping force necessary to the 
suspension to mitigate any kind of unwanted vibration which oc- 
curs in many systems such as seismic, wind induced frequency, 
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prosthesis control, washing machine, vehicle suspension. When- 
ever a problem happens in the system, the magneto-rheological 
damper has the ability to react as a passive damping and lithium 
stearate. These three components of an MR fluid ensure its mag- 
neto rheological characteristics. Modifying any one component 
will lead the variation of rheological and magneto rheological 
properties of the MR fluid. An important combination of all the 
three elements is important to reach the required properties of MR 
fluid. 


2.3 Composition of MR fluid 


The magneto rheological fluid is composed by three elements, 
base fluid, metal particle and additives. 

Base fluid is a passive or non-magnetic carrier fluid in which the 
micro-sized iron particles are suspended. The base fluid must pos- 
sess natural lubrication and damping specifications. In design and 
good application of MR fluid; base fluid should have a low vis- 
cosity and not change with temperature. It is very important that 
variations of viscosity caused by magnetic field become high 
compared to variation caused by natural physical parameters. Ac- 
cording to the presence of suspended particles base fluid becomes 
dense. Famous utilized base fluids are hydrocarbon oils, mineral 
oils and Silicon oils[24]. 

Metal particle for the best used of this technology a type of 
metal particles which can provides a fast and easy response of 
magnetization is needed. Metal particles used in the MR- technol- 
ogy are very small. Size of the particle is around range of | micro- 
meter (tum) to 7m. The most used and available on market metal 
particles are carbonyl iron, powder iron and iron cobalt alloys. The 
properties of these materials provide high magnetic saturation ac- 
cording to that they are is capable for, a good magnetizing chain 
structure. The concentration of magnetic particles in base fluid 
may reach to 50%. 

Additives to have a better control of MR fluid properties some 
amounts of additives are required. It has the role of stabilization and con- 
trol the viscosity. Additives are composed by stabilizers and surfactant; 
surfactants ensure the decrement rate of settling of the metal particles. 
While the functions of additives are to control and monitoring the vis- 
cosity of the fluid, keep friction between the iron metal particles and to 
minimize the rate of reduction of the fluid due to long period use of the 
fluid consequently additives avoid wear of fluid and increase it service 
life[25]. Most utilized additives are ferrous oleate 


3. Magneto rheological dampers 


Based on specific application, configuration, desired load 
MR damper have been designed and developed in different sizes. 
The configuration of a damper can be classified into three main 
categories based on operation mode, hydraulic housing (cylinder), 
piston (or MR control valve) structure. Damper is composed by 
cylinder housing filled with the MR fluid and a control valve (pis- 
ton) which generates damping force while it is in motion. Gener- 
ally the control valve is embedded with solenoid which is the mag- 
netic circuit. According to the operating principle modes MR 
dampers are categorized into four fundamental types: flow-mode 
dampers, shear-mode dampers, squeeze-mode dampers and 
pinch-mode dampers. Dampers which are using more than one 


MR fluid operation mode are called hybrid or mixed mode devices. 
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Every type has brought particular interest, and it might use various 
objectives. 


3.1 Mono tube damper 


Mono-tube dampers represent the commonly used magneto 
rheological devices. Due to its simplicity design and less internal 
parts mono tube damper remain the best choice for different ap- 
plication. It is widely used in passenger vehicle suspensions. 
Based on the standpoint of design and implementation of MR 
dampers, vehicle suspension using MR dampers, mono-tube pre- 
sents a significant benefit of no having electromechanical valve 
and minimum number of moving elements. The working principle, 
the design setting of a characteristically gas-charged MR damper 
is the same as that of a natural valve-based damper. The floating 
piston serves on separation of MR fluid from high pressure gas 
chamber. The role of main piston is to separate the quantity level 
of magneto rheological fluid into compression chamber and re- 
bound the fluid volume between piston road and main assembly. 
As it shown in Figure 6 the piston assembly embeds an annular 
gap (annulus) to allow the fluid to move between the chambers 
while in motion. The relative movement of the wheel (unsprung 
mass) and the body (sprung mass) conduct the fluid flow between 
the chambers through the annular gap in the piston. The big pres- 
sure gas in the chamber under the floating piston is needed for 
cavitation-free operation, but, it has an opponent effect on cylin- 
der inner surface complete necessaries, sealing and friction[3]. 


MR Fluid Reservoir 
Piston Guide 
Piston 


Accumulator Piston 


Compressed Gas Reservoir 


Figure 6. Mono tube MR damper 
3.2 Twin-tube Damper 


As it called the twin-tube MR damper is one of MR damper 
that possess two fluid reservoirs, one located inner the other one 
as it is indicated on the Figure7. In this type of setting, the damper 
possesses two housing an inner and outer housing[26]. The inner 
housing ensures the guidance of the piston rod assembly, in ex- 
actly the similar method as in a mono-tube damper. The volume 
stayed by the inner housing is assigned to as the inner reservoir. 
Equally, the volume that is determined by the space between the 
inner housing and the outer housing is assigned to as the outer 
reservoir. Hence no air pocket exit because inner reservoir is 
filled by magneto-rheological fluid[3]. 


Piston Rod Foot Valve Assembly 


Inner Housing Outer Housing 


Figure 7 Twin tube MR damper 
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To arrange variations in volume caused by piston motion, the use 
of outer reservoir partially filled is necessary. For, the outer hous- 
ing in a twin-tube damper shows the similar objective as the pneu- 
matic accumulator system in mono-tube dampers. Practically, the 
valve assembly named foot valve is mounted to the bottom of the 
inner housing to control the flow of fluid between the two tubes. 
Whenever the piston rod, move into the damper MR fluid move 
from the inner tube to the outer tube through the compression of 
foot valve assembly. The quantity of fluid that flows from the in- 
ner reservoir into the outer reservoir is the same as the volume 
moved by the piston rod as it enters the inner housing. When the 
piston rod is returned from the damper, MR fluid flows from the 
outer reservoir into the inner housing tube through the return valve, 
which is also part of the foot valve assembly part. 

system with known operational parameters depending on specifi- 
cations of the magneto rheological fluid when a magnetic field is 
not exposed on it[21]. 


3.3 Double ended tube MR damper 


The last type of MR damper is named a double-ended 
damper because a piston rod of equal diameter exceeds from both 
ends of the damper housing. Figure8 indicates a cross section view 
of characteristically double ended MR damper. Because of no var- 
iations in volume exists; when the piston rod moves relatively to 
the damper body, in the other word the double-ended damper does 
not need an accumulator system. This kind of MR damper is ap- 
plied in gun coil system, bicycles, seismic vibration avoidance, 
and wind drift[21]. 


MR Fluid Reservoir 


Rear Piston Rod 


Front Piston Rod 


Piston Approximate Flux Path 


Figure 8 Double-end tube MR damper 


4. MR damper fundamental models 


The fundamental model of MR fluid is based on Bingham 
solid, which possess the variable yield strength (elasticity)[16]. In 
Bingham model, fluid flow is provided by Bingham’s mathemat- 
ical equations, 


T= Ty(H) + ny (2) 


T<Ty 


tT represents the fluid stress, t, denoted the field dependent 
yield stress, H defines the magnetic field, y is defined as the fluid 
shear rate, and 7 denotes the plastic viscosity; it means, the vis- 
cosity if H=0. Under the fluid’s yield stress (pre-yield state), the 
fluid displays visco elastic characteristics[27]. 


43 


4.1 Flow mode 


The working principle of this mode is based on fluid which moves 
through the two fixed surfaces and magnetic field is exposed per- 
pendicularly to the direction of flow as shown in Figure 9. By ad- 
justing the intensity of magnetic field, the viscosity of fluid is con- 
trolled. This mode of MR fluid technology is utilized in different 
categories of dampers and vibration absorbers and has wide appli- 
cation in automotive industry. In classic dampers, fluid moves 
from the inner reservoir (cylinder) to outer reservoir (cylinder) 
through the foot valves. The air and fluid is located in outer. When 
piston moves down, the fluid volume level in outer cylinder raises 
and free air behaves as a compressed medium and generates 
damping force. When the piston spread out, direction of flow gets 
reverse MR dampers are slightly opponent from the tradition 
dampers. This type technology does not have the valve like clas- 
sical dampers. The most application of this mode is the suspension 
of highway vehicles. In this kind of damper there is an annular 
orifice passage, through which MR fluid can be transferred from 
one chamber to another. In normal operation mode ( in absence of 
magnetic field)condition, fluid can flow through the orifice but 
when suspension is needed, the solenoid is energized and current 
starts to pass through the coils windings, and generates the inten- 
sity of magnetic field necessary to provide a damping force[24]. 


South(-) 


pressure 


tg FowQ) = 


| North(+) 
Applied field 


Figure 9 Flow mode 


As it is shown in Figure 9 above, pressure variation in this 
mode is the calculated as sum of pressure generated due to fluid 
viscosity AP and pressure developed due to the magnetic field 
AP,. 
AP = AP, + AP, = 12nQL/(g?w) + ctyL/g (3) 

Where AP, denoted viscous component of pressure drop, 
AP, field dependant of induced yield stress part of pressure drop, 
n plastic viscosity (H=0), g, fluid gap, W defines Width of fluid 
flow orifice, Q pressure fluid, t, defines Field dependant yield 
stress found in Magneto rheological fluid, c, constant, L length of 
fluid flow orifice. 


4.2 Shear mode 


The operation of this mode is based on the fluid which moves be- 
tween surfaces having same motion and a magnetic field is sub- 
jected perpendicular to the direction of flow as shown in Figure 
10 [28]. Shear mode of MR fluid technology is utilized in different 
type of brakes and clutches of automotive. In this mode, the total 
shear force generated is a summation of the force generated based 
on viscosity of the fluid (F,) and the force generated based on the 
magnetic field (F,,). 
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Figure 10 shear mode 


The design of squeeze mode Figurel1 is simple and pos- 
sesses the minimum moving elements hence, less maintenance is 
desired, and service life is increased. The main disadvantage of 
MR fluid technology is that the MR fluid becomes dense after 
elongated use and needs to be changed, also the weight of MR 
fluid is high caused based the high density of iron metal particles. 
Presently automotive industries are applying this system. The ap- 
plication of this technology is wide such as, medical engineering, 
aerospace and polishing. There is a large range of research in MR 
fluid technology. MRF systems have to be designed more delicate, 
probably by introducing the uses of sensors and feedback technol- 
ogy i.e. closed loop systems. The life range of MR fluid devices 
in terms of the total energy dissipated from the equipment need to 
be increased. An improved MRF technology would make it the 
smart technology of future[28]. 


displacement force 


Figure 11 Squeeze mode 


5. Mechanical model of MR damper 


MR dampers are called semi-active elements of suspension 
systems whose damping specifications may be changed in real 
time based on the ability to change the viscosity of MR fluid as 
well its resistance to circulate by using of magnetic field. The ben- 
efit of MR fluid is that can permit the changes in the damping 
force by varying the current subjected to the damper. Thereby the 
hysteretic model characteristics of MR damper are based on cur- 
rent dependent and a function of excitation amplitude. It is im- 
portant to analyze and model the MR damper with appropriate 
methods in order to predict it behavior in various magnetic fields 
or excitations. This approach becomes important aspect because 
before the production of a MR damper it is significant to design 
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and choose the right parameters that will reach the behavior of the 
damper[29]. 


5.1 Bingham model 


This Bingham model investigates and study plastic charac- 
teristic and the behavior of magneto-rheological and electro rheo- 
logical fluids. Bingham model ensure that a body remain as solid 
till a little yield stress is surpassed then after it shows a linear re- 
lationship between the stress and the rate of shear deformation[30]. 
Some researchers like Stanway et al. investigated a mechanical 
model based on the Bingham plastic model to exhibits the ER 
damping system. As it indicated in Figure 12 the model is based 
on a combination between coulomb friction and viscous dashpot. 
and is always utilized to analyze the characteristic or behavior of 
MR dampers[31]. 
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Figure 12 Bingham model 


The force provided by mechanical element setting is coming 
from equilibrium concepts; it is expressed in the equation below: 
F(t) = cox + fosgn(X) + fo (5) 

Where, co denoted as constant of damping coefficient, x de- 
fines the velocity of the external excitation fc denoted the fric- 
tional force and fo is the force offset related with the availability 
of an accumulator considering that this element has linear behav- 
ior and low stiffness. To investigate the flow of fluid, this mathe- 
matic expression model can be derived. 


5.2 Bouc-Wen model 


The hysteresis Bouc-Wen equation is among different mod- 
els and the most used to study and simulate various method of 
expanding the damping system by changing the hysteresis. It has 
been invented by Bouc in 1971[32] and later developed by Wen 
in 1976[33] who has showed different methodologies of this 
model to show a wide kind of hysteretic patterns. According to 
this interesting behavior, the model was utilized to analyze differ- 
ent nonlinear hysteretic systems like MR damper and hysteretic 
isolators. However, MR dampers show a good typical behavior, 
the final hysteretic setting relies on certain particular specification 
features related to accumulator, damper geometry and so, on. So, 
the model should be suitable to have the real MR damper charac- 
teristic and several changes of the Bouc-wen model were designed 
and developed to simulate properly MR damper. Based on Bouc- 
Wen model different procedure to study the simulation the MR 
damper has been proposed; these models are: Classic or simple 
Bouc-Wen model; 
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Advanced Bouc—Wen model or Spencer model; 
@) Shear model Bouc-Wen model 

(3) Bouc—Wen model for wide scale MR dampers; 
(4) Inverse Bouc—Wen model which is current 
dependant; 

(S) Current-frequency—amplitude-dependent; 

(6) Asymmetrical Bouc-Wen model. 

As, it is indicated in the Figurel3, the simple Bouc-Wen 
model is composed by three elements: a spring, a Bouc-Wen block 
and dashpot, mounted a parallel setting. In my project the simple 
model of Bouc-Wen will be used to simulate the MR damper be- 
haviors. In this figure nonlinearity of the system is situated in the 
Bouc-Wen block, which is able to track the behavior of MR 
damper: 
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Figure 13 Bouc-wen model 


The model investigates different characteristics of MR 
damper by comparing the performance of this model with other 
parametric models. Based to mechanical setting indicated in Fig- 
ure13, the equation below express the damping force. 


F(t) =cox+ k(x —xX9) + az (6) 


Where co denoted the constant of viscous coefficient, k de- 
noted the constant of stiffness coefficient and z is an evolutionary 
related variable with Bouc-Wen block and it ruled by: 


z= —ylx|z|z|""+ — Belz" + Ax (7) 

Xo is denoted the displacement at initial state, which permit 
the presence accumulator in the mechanism system. Constant co, 
k, a, 6, y, A and n are generally and the Bouc-Wen shape parame- 
ters or behaviors. They are the functions of the current subjected 
to the MR damper, frequency of vibration and amplitude. Nonlin- 
ear shape hysteretic curve can be controlled by varying the 
amounts of the Bouc-Wen block parameters permitting to monitor 
the smoothness in the state of transition between pre yield and post 
yield. Based to the ability to take the hysteretic element of the MR 
damper characteristics, the classic Bouc-Wen model is consid- 
ered a best rational model of mathematical and numerical simula- 
tion compare to Bingham model. It is better in analysis of hysteric 
behavior. The use of Bouc-Wen hysteric model, lead on a good 
simulation of MR damper, test and implementation of its par- 
ties[35]. 
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6. Semi-active vehicle suspension control strategies 


Semi-active suspension is considered as the union of active 
and passive suspension. Since its invention, the damping force 
was adjusted by the variation of orifice area which leads on the 
change of flow resistance. MR damper is a kind of semi-active 
damper where the motion of MR fluid is adjusted by changing the 
value of current applied and thereby varies the level of damping 
force. As other control strategies of semi-active suspension damp- 
ers the MR damper is capable to provide a damping system even 
if the current become zero. The damper is able of generating huge 
control forces by changing behaviors like stiffness and damping 
coefficient through the variation of magnetic field and thus con- 
trols the response of large scaled structures like big buildings [36]. 
Control strategies of this system are to give a good stability han- 
dling and a smooth riding. Many studies and researches have been 
conducted to optimize the trade-off between road handling and 
ride comfort. These control strategies are categorized into tech- 
niques based on, nonlinear, linear and intelligent control method- 
ologies. 


6.1 Skyhook control 


Skyhook control strategy is among the semi-active control strate- 
gies and the most commonly utilized that changes the damping 
force to imitate the damping effect of a passive damper in its op- 
timal setting, mounted between the fixed reference and sprung 
mass (vehicle body) as showed in Figure 14 [58]. 
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Figure14 skyhook control strategy 


It is characteristically categorized as continuous control 
skyhook and on-off skyhook system. On-off Skyhook controller 
is always easier designed and well appropriated for the industrial 
era [37]. To implement the control strategies of On-off law equa- 
tion below should be governed. 


if V.Vz =0 
otherwise 


Gg] (8) 


Cnin 

Where, V2 is defined the absolute velocity of the body, and 

V, denoted as the relative velocity of the sprung mass (body) and 
unsprung mass (wheel) across the suspension. It is explained that 
if relative velocity of sprung mass (body) and the wheel have the 
similar direction, Hence, the maximum damping force must be 
subjected to mitigate the body acceleration. Meanwhile, if the two 
velocities (sprung velocity and unsprung velocity) are in the odd 
directions, the damping force must be set at low level to mitigate 
body acceleration. To achieve the skyhook controller strategy, ab- 
solute velocity of the body should be measured. Skyhook strategy 
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has been proven the optimization of ride behavior by minimizing 
vertical body acceleration. Skyhook control is good approach to min- 
imize the vertical vehicle dynamic. But, according to its emphasis 
on optimizing a smooth ride, the road handling performance is 
very small. To prevent this issue ground-hook control strategy has 
been designed and proposed to optimize road handling abilities 
for controllable suspensions. The theory is based on the virtual 
damper which mounted on a fixed ground instead of sky which 
will give consideration to road holding ability to suspension 
mechanism[38].By comparing the skyhook semi active controller 
with passive suspension; models are simulated using a square 
bump as road profile, results show that the skyhook presents a 
great contribution to comfort of vehicle occupants as it is shown 
to the figure. 15 Where minimum body acceleration is indicating 
ride comfort 


0.2 T 
— Passive suspension 
— Skyhook controller 


Body acceleration (m/s2) 


0.15 l 1 Ll Ll | 
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Fig: 15. Body acceleration comparison skyhook vs Passive 


6.2. Ground hook control 


As it is indicated in the Figure 16, the ground hook strategy 
is opposite from the skyhook model. In this model, the damper is 
mounted to the wheel rather than the body. Under the ground hook 
setting, the focus changes from the sprung-mass to the unsprung- 


mass [60]. 
ae 
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Figure 16.Ground hook strategies 


As it indicated in this Figurel6, in this setting the focus 
shifts from the sprung mass to unsprung mass. While skyhook 
work well in controlling the vertical body mass, ground hook is 
good in controlling the unsprung mass frequencies. The idea of 
ground hook is the same as the skyhook strategy, the difference is 
that it is ground hook control the wheel vibration [39]. 


Faq = ‘gnd V; 


gee >0 
Fog = 0 


VV. <0 (9) 
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6.3. Sliding control mode 


Sliding control mode is a strong control strategy with various sig- 
nificant functionalities like insensitivity to linked behavior 
changes and disturbance isolation. Depending on its applications, 
sliding mode controller design showed a great potential in keeping 
the stability and optimization of in modeling and simulation the 
uncertainties. Sliding controller is constituted by of a nominal part, 
the same closed loop theorem, and additional conditions objected 
to face with uncertainty model. Reference model of Sliding Model 
Control is consist on forming a ‘sliding surface’ to track errors, 
the tracking should be equivalent all time[40]. 

Sliding-mode design has various functionalities that make it 
an advantageous option to remove nonlinearity tracking process, 
however; it presents a great importance in designing and modeling 
of uncertainties systems. It can be explained in many ways[41]. 
The first way is to consider that the system to process has two step 
procedures. Primo is the sliding surface which the desired state of 
the system. Secondly, a control system to conduct toward a sliding 
surface in the other word a system to track the sliding surface sys- 
tem, then after switching on that surface. In theory sliding mode 
will effectively approach the control purpose, the cause of phase 
surface[42]. The control action design can be pretended according 
on various strategies, a straightforward one which is based on de- 
termining a Lyapunov-like condition that make the sliding surface 
to be an attractive region for the state vector trajectories[43]. 

Reference: the idea of reference is based on development of 
control algorithm necessary to force the plant to virtual dynamic 
model. So, to consider a reference model is very useful. In model 
reference control for instance a set point regulation problem is 
changed into a tracking problem by immerging a reference model 
to filter the given configuration point value and to provide output 
which is variable with time as virtual response of tracking control 
mechanism[42]. 

Based to the sliding mode control (SMC) regulation, as it 
is indicated on the Figure17a and b, an asymptotically fixed slid- 
ing mode state can be derived from dynamic errors between the 
MR quarter suspension of vehicle dynamic system and the trans- 
formed skyhook reference model approach[44]. 
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Figure 17. Sliding control mode 


According to skyhook reference model wheel (unsprung 
mass), the displacement x, of MR quarter-vehicle suspension can 
be directly utilized as road variation input of the skyhook refer- 
ence model system, because the real tire stiffness is typically 5 
times of the suspension spring stiffness. Based on hypothesis 
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above the dynamic equation of skyhook reference is expressed 
as[45]: 


MsoXso a Co(Xs0 = Xy) + CsoXso + ks (Xs0 = Mu) =0 (10) 


Where m, 9 denoted as reference model sprung mass, 
which is equivalent to unloaded sprung mass of the semi-active 
quarter-vehicle suspension, Cs 9 denotes skyhook damping con- 
stant, and X,q are defined as body mass displacement of refer- 
ence model, k, and cg defined as spring stiffness of the passive 
quarter-vehicle suspension and damping constant coefficient, re- 
spectively. The sliding mode s of MR quarter-vehicle system is 
expressed as this equation below: 


s = é,(t) + Ae, (t) (11) 


where A>0 defines the convergence rate of sliding mode sur- 
face, e,(t)=x,(t)-x,9(t) defines the errors between body mass 
displacement of MR suspension quarter car dynamic and dynamic 


system of skyhook reference and e(t) become 0 if t tends to infinite. 


Moreover, Lyapunov function is expressed as V=s?/2. To make 
sure that the system state asymptotically adjusts to the sliding sur- 
face mode, it is denoted that 


V =ss < —4Is| (12) 


Where @, is denoted as positive constant. If the effect of 
gravity is neglected, the equation (13), is gotten 


Mk Khe — oy) + Fy = mg = 0 


fs 13 

tes +g kee =—%,) — Fy yg =0 a) 
We get 

X, = [—k,(%5 — x,) — Fal/ms (14) 


By combining Eq. (13) and Eq. (14), equation (15) is ob- 


tained 
s= [—k5(Xs =i Fql/m,; — Xoq + Ae, (t) (15) 


Base on the stability functionality of sliding mode control 
strategy, $s=0, the virtual control damping force FdO of MR 


quarter car suspension can be derived from Eq. (15) 
Fao = —ks (Xs — Xu) — MsoXso + MsoAes (16) 


To guarantee that yielded damping force Fg of the MR 
damper can quickly track its virtual control damping forces Fo, 
by using switch control law, the damping control force is denoted 
F, = Fag — K + sgn(s) (17) 

Where K defines the constant gain of SMC which will be 
derived later and sgn(s) is limit function which vary between - 
1 and 1, If s<O and s>0, respectively. Based on V= Ss, F, to re- 
place Fy in Eq. (15) integrate it into V= $s and therefore, obtain 
Vis gotten 
V= (-k,x5/M, + k.xy/M, — F,/Ms — Xsq + A€s)S 

Combining Eq. (16) and (18) V, is obtained 
V= [/m, — 1/m59)(—k x5 + ksXy — Fao) + K + sgn(s)/ms] (19) 


(18) 


Based on the Lyapunov stability of Sliding control mode, 
Eq. (19) is integrated into Eq. (14) and obtain 
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1 
= =) (-k x5 + KX, -_ Fao) 
+K -sgn(s)/ms 


s<-als| 


(20) 


[Ge 


As considering s =s: sgn(s), Hence Eq. (20) is governed as 


[br = oa) Chats tise Fao) 5 2 age seas) (21) 
+K + sgn(s)/m, 
for s>0, Eq. (21) can be simplified as 

K < —mg@ + (1 — m,/mg9) (Kk 5X5 — KsXy + Fao) (22) 
for s<0, Eq. (21) can be simplified as 

K < —m,¢ + (1 — m,/Ms59) (kg X5 + ksXy — Fao) (23) 


If the real load of vehicle is considered as usually variable 7 
is expressed as the change rate of vehicle body mass (sprung 
mass) m, to its unloaded sprung mass Mg, in the other word 
S=M,/Ms9=1. By using Eq. (22) and (23), K of sliding control 
mode is obtained 


K= —M5oND — (n ~ 1) (lFaol + ks |x| + Kglxy1) (24) 


Therefore, the saturation function saf(s) is used to replace 
the sign function sgn(s) in Eq. (20) for limiting chattering fre- 
quency which may happen in the system, and takes s/e and 
sgn(s) when s <eand s > ¢, accordingly. Eq. (20) can be shown 
as 


F, = Fag — Ksat(s) (25) 
Where ¢ denotes positive constant. Fig18, indicates the ride 
comfort between passive suspension and semi active suspension 


in sliding mode control strategy. The improvement is more than 
35% compare to passive suspension. 
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Figurel8. Body acceleration comparison of sliding mode and passive 


3.4 Fuzzy logic control strategy 


Fuzzy logic control strategy is a type of many values of 
logic. In many applications of control of nonlinear systems, this 
approach developed by mathematicians has been appreciated and 
successful used. Since the last decades the implementation of 
fuzzy logic has been growth due to the various researches which 
showed its potential on controlling the nonlinearities systems, 


Musabyimana et al. 


Engineering Perspective 1 (2): 38-51, 2021 


Fuzzy logic has been inverted by Zadeh [46] in 1965 , ten year 
after its invention , Mamdani [47] developed and showed the util- 
ity of Fuzzy Logic Control (FLC),in modeling and control of 
small steam engine’s suspension to replace the tradition suspen- 
sion of damper and spring of vehicle by active suspension. The 
result of Fuzzy logic model demonstrated a big improvement of 
riding comfort and road handling in various conditions [1]. Pres- 
ently fuzzy logical is applied in many areas such as dynamic pro- 
cess, washing machine and in automation industries. The theory 
of fuzzy control is based on operations rules that are human can 
read and represent human’s knowledge about how to control and 
monitoring the process. The rules are defined in terms of fuzzy 
involvement. It is composed by three main elements namely input, 
process output and a reference signal. For the suspension vertical 
dynamics application, input is the suspension damping control and 
process output is a rate of velocities between sprung and unsprung 
masses [49]. 

Generally, input is supplied to Fuzzifier which changes 
crisp data input into fuzzified data which is then analyzed pro- 
cessed in Fuzzy inference function block using fuzzy rules. This 
fuzzy output is then converted to crisp output using Defuzzfier. 
Fig 19 indicates schematic of Fuzzy control methodology which is 
the flowchart of control strategy. Fuzzy logic approach control has 
many features that make it attractive for semi-active suspension. 
There are a many nonlinear devices like tires shock absorbers 
springs and many rubber components. All these non-linearities 
can be made linear having significant impact on the precision of 
model. 


Fuzzy rules 


Fuzzy 
inference 
function 


Input Output 


= ==}| Defuzzifier 


Figure 19 fuzzy control methodology 


Fuzzy logic algorithm controls many loops variation in var- 
ious non-linear systems. it has two inputs, namely suspension de- 
flection and suspension velocity. the controller calculates those 
inputs and generate the desired damping force based on inference 
rules subjected on it. These control systems embed the human as 
thinking through the application of fuzzy configurations and lin- 
guistic variables related by a configuration of IF-THEN fuzzy 
tules. The results show figure 20 and 21 that fuzzy logic can dra- 
matically reduce the tradeoff between road handling and ride com- 
fort. However, it can’t completely mitigate this odd. 
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Figure 20 Body acceleration comparison of semi active and passive 
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Figure21 Suspension deflection comparison of semi active and passive 


3.5 Linear Quadratic Gaussian 


Linear-Quadratic-Gaussian (LQG) control is a robust strat- 
egy which is based on an advanced state-space algorithm to design 
optimal dynamic regulators. It is based on noise measurement of 
the process and requires a state-space model of the plant[50]. This 
control algorithm was less used in control of semi-active vehicle 
suspension based on MR dampers, thus it was largely applied in 
control and monitoring of active suspension to mitigate the vibra- 
tions. To examine the effectiveness of semi active suspension con- 
troller, the LQG control strategy is utilized to calculate the desired 
damping force fy necessary to generate a smooth ride and comfort. 
LQG is composed by Kalman filter gain Kyand of an optimal state- 
feedback, gain K named Linear-Quadratic-Regulator (LQR). The 
optimal control has been designed to find a controller which may 
perform better. The optimal controller LQR for the vehicle sus- 
pension control is the most utilized in background of various re- 
searches of vehicle suspension control. The strength of LQR 
method is that in using it the factors of the optimization index can 
be weighted based to the designer’s wish or other constraints. The 
LQR approach is utilized to optimize the road holding and the 
drive comfort for vehicle suspension. LQR control approach re- 
main better in control of linear active suspension mechanism [50]. 

In modeling and control of linear active suspension, the 
equation below denoted the system 


x(t) = Ax(t) + Bu(t) + f@) (26) 


A state variable of feedback regulator of the system is con- 
sidered as 
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u(t) = —Kx(t) (27) 


K denotes the state feedback constant gain matrix (control- 
ler). Hence, a new system is obtained as 


x(t) = (A— BK)X(E) + FO) (28) 


The perform procedure serves on defining input of control 
U, which reduce the performance index. J is performance index 
show the optimization required behavior as well as the controller 
input mitigation. The LQR optimal controller of provided system 
is defined as controller design which reduces the performance in- 
dex below. 


J= fo XOEX) +u(t)’Ru(t)) dt (29) 


Where u defines state input, R and Q denoted the positive 
weighting matrices. 


Controller gain K is defined by: 
k= “HP (30) 


Where P, is denoted as the matrix which is studied as solu- 
tion of the Algebraic Riccati Equation (ARE) as follow as 
PA+A™P — PBR '1B7P+Q=0 G31) 

To evaluate the Kalman filter gain Kf, the solution of the 
following Algebric Riccati Equation is obtained as S 


AS + SAT + GWGT — SC'V-105 +Q =0 (32) 
Then, calculating Ky like that, 

Ky =SC'y"1 (33) 
The Kalman filter for the LQR controller becomes 

R(t) = AR(t) + Bu(t) + K-[y(@) — CR(t) — Du(t)] (34) 
The equation below is considered as system 

x=Ax+Bf,+Gw (35) 


Where w(f) is defines as a fixed Gaussian white noise pro- 
cess. Then, system (29) is denotes as linear system applied to a 
Gaussian noise. By solving stochastic equation of optimal control 
issue, and measuring the feedback system the equation (36) is got- 
ten. 


Y = Cx(t) + Du(t) + v(t) (36) 


Where v(t) defines a random signal which defines the meas- 
urement noise. The plant and measurement noise are remained 
both white noise and simultaneous. As it is shown below, 


Efww}}=W,  Efvv'}=V, E{wv'}=0 (37) 


The LQR gain K defined above are utilized with a Kalman 
filter gain Kyto get state estimator. It was determined that the body 
acceleration x» and the wheel acceleration x,, measurements are 
available in addition to the suspension working space x. In other 
words, the output vector is assumed to be Y = [Xp - Xw, Xp, Xv]. The 
measured parameters are utilized to optimize LQG controller of 
suspension because the cost accelerometers are relatively low, and 
LQ was not robust enough to perform a better control. 


4. Conclusion 
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Semi active control strategies based on Magneto-rheologi- 
cal damper are studied and reviewed in details, through the Bouc- 
wen model for MR damper, the simulation was conducted based 
on quarter car model passing on trapezoidal bump, the control im- 
provements are just subjective and controllers should act with a 
compromise between odds behaviors which is ride comfort and 
vehicle stability. Simulation results for various control strategies 
like skyhook, groundhook, slide mode controller, fuzzy logic and 
Linear quadratic Gaussian for Bouc-Wen magne-torheological 
damper, showed that the semi active suspension acts well on the 
reduction of compromise between comfort and stability, however 
it can’t mitigate this compromise at 100% this odd. The authors 
suggest that combining more than two control strategies could 
dramatically reduce this compromise issue between comfort and 
vehicle stability. For the future scope the use of artificial neuro- 
networks in control strategy could bring lower odd between ride 
comfort and road handling. 
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